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Abstract. Entanglement dynamics for a couple of two-level atoms interacting with 
independent structured reservoirs is studied using a non-perturbative approach. It is shown 
that the revival of atom entanglement is not necessarily accompanied by the sudden death 
of reservoir entanglement, and vice versa. In fact, atom entanglement can revive before, 
simultaneously or even after the disentanglement of reservoirs. Using a novel method based 
on the population analysis for the excited atomic state, we present the quantitative criteria 
for the revival and death phenomena. For giving a more physically intuitive insight, the 
quasimode Hamiltonian method is applied. Our quantitative analysis is helpful for the practical 
engineering of entanglement. 
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1. Introduction 

Entanglement is highly relevant to the fundamental issues of quantum mechanics and plays 
a central role in the application of quantum information [[TJ. In recent years, it has been 
observed that two qubits interacting with independent reservoirs experience disentanglement 
in a finite time in spite of the asymptotical decoherence. The corresponding investigations 
for this phenomenon, called entanglement sudden death (ESD), give notable results both 
theoretically [|2l [31 H |5l [6l 13 and experimentally (HI 13. Nevertheless, most of previous 
works concentrate on the dynamical evolution of a bipartite system, while the understanding 
of how the information is transferred is to be explored in depth, especially for non-Markovian 
environment. 

Recently, people begin to analyze the dynamics of an extended system which 
incorporates two qubits and their respective reservoirs and successfully explain the process 
of information transmission in Markovian regime ifTOl . Thanks to these significant 
improvements and using the non-perturbative method [fTTl [T2l [13)1 . we extend the study of 
this composite system into non-Markovian regime where the reservoir presents some internal 
structure. 

First of all, we study the entanglement dynamics of two atoms and two reservoirs. 
Expectedly, when the atom entanglement is depleted permanently, it is completely transferred 
to the reservoir entanglement irrespective of the non-Markovian memory effect. The 
reservoir entanglement can not only exhibit entanglement sudden birth (ESB) but also ESD. 
Surprisingly, the revival of atom entanglement [ fT4| is not always accompanied by the 
disentanglement of reservoirs, and vice versa. Entanglement of atom pairs can revive before, 
simultaneously or even after the ESD of reservoirs and these phenomena are independent of 
the relative strength of atom-reservoir interaction. 

In addition, we analyze the effective bipartite entanglement and multipartite 
entanglement within the system and thus give a comprehensive interpretation of the 
information transmission in non-Markovian regime. Applying the population analysis for 
the excited atomic state, we present quantitative criteria, which have been left obscure for 
quite a long time, of the occurrence and exact numbers for revivals of atom entanglement 
and death of reservoir entanglement. Our analytic and numerical analysis enables one to 
precisely manipulate the entanglement based on the atoms or/and quantum dots in high Q 
cavities [|26ll27ll28ll29ll and thus have potential importance in application. 

Finally, by transforming the true mode Hamiltonian [|T2l [T9ll to the quasimode form [[T2| 
which explicitly separates the memory effect and the damping effect of the reservoir, we give 
a more physically intuitive insight into the non-Markovian phenomena. 

2. Theoretical framework 

Under the rotating wave approximation, the truemode Hamiltonian of the single-body system 
(atom and its reservoir) is (/i = 1) [|20ll2T1l 



H = Ho + Hi 



(1) 
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da>k(Okb\oJk)b{oJk), (2a) 

CO 

Xco 
da)tg(oJk)cr+b(a)k) + H.c. (2b) 
CO 

Here, cr+ and coq are the inversion operators and transition frequency of the atom, and 
b(oJk), b^(oJk) are the annihilation and creation operators of the field mode of the reservoir 
with the eigenfrequency cok- 

For an initial state of the form \e) (S) 10)^ with 10)^ = Y\k \(^k)r^ the time evolution of the 
single-body system is 



Xco 
da;,c^,(Oe-''^>'' \g)\h)r, 
CO 



(3) 



where \\k)r is the state of the reservoir with only one exciton in the kth mode. According to 
the Schrodinger equation, the equations for the probability amplitudes take the forms 



ico(t)= r c,,(t)g(a;k)e-'("^-"'')'da;k, (4) 



ic,,(t) = g*(a;k)e'('^^-'^'')'co(t). (5) 

Eliminating the coefficients c^j(0 by integrating ^ and substituting the result into (31), one 
yields 

co(t) = - [ co(ti)f(t - ti)dtu (6) 
Jo 

where the correlation function takes the form 

Xoo 
da;kJ(c^k)e-'^"^-'^°^^'-"\ (7) 
CO 

Suppose the atom interacting resonantly with the reservoir with Lorentzian spectral density 
J(cok) = \g(cok)\^ = W^A/n[(cok - c^o)' + A\ (8) 

by employing Fourier transform and residue theorem, we get the explicit form f(t - ti) = 
W^exp(-A \t- ti\), where W is the transition strength and the quantity l/A is the reservoir 
correlation time. 

To construct a concrete physical insight, we adopt a nonperturbative method called 
pseudomode approach [fTTl ITU [T3]| in the following analysis. Equation ^ can be restated 
as 

co(t) = -iWb(t), (9) 
b(t) = -Ab(t) - iWco(t), (10) 



where 



b(t) 



-iW r Co(ti)e*'Mti 
Jo 
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is the pseudomode amplitude. Typically, there are two regimes lfT3 : weak-coupling regime 
{A > 2W), where the behavior of the single-body system is Markovian and irreversible 
decay occurs, and strong-coupling regime {A < 2W), where non-Markovian dynamics occurs 
accompanied by an oscillatory reversible decay and a structured rather than a flat reservoir 
situation applies. We will limit our considerations to the latter case. 
Define the normalized collective excited state of the reservoir as 

Xoo 
c.,a)/c2(0e"'"^'IU),da>,. (11) 
CO 

one can rewrite dS]) as 

1^(0) = ci(Ok) 10), + C2(0 1^)11)., (12) 

where 



ci{t) = e 



-(iwo+/}/2)t 



dt\ A . Idt 



with d = {AW^ - A^yi^. We note that ci(0 = Co(0 exp(-ia;ot), C2{t) = (1 - \c^{t)\y^ which 
can be calculated directly from the definition (fTTI) and 

Z7(0 = -2i^e-'^'/2^in(^ 
d \ 2 

The atom and its reservoir now evolve as an effective two-qubit system. 

Now, we study the joint evolution of two identical single-body systems initially in the 
global state 

|0o) = {a \g), \g)2 +/3\e), \e)2) |0),, 10),, , (13) 

where the real non-negative parameters a and /3 satisfy +/3^ = 1 and / (/ = 1, 2) denotes the 
zth single-body system. The evolution of the composite system reads 

10(0) =a\g),\g)2\0),^\0\^+/3\ipm\n(t)), (14) 

where \(pi(t))(i = 1,2) represents the single-body evolution and can be determined by (fT2l) . 
Employing the density matrix p{t) = |0(O) (0(01 and the definition of concurrence 1221, we 
can write down the concurrence of entanglement for different partitions: a\® a2, r^® r2 and 
a\ (8> ri, respectively, 

Ca,cn{t) = max{0,2;e|cip(a-y8|c2h}, (15a) 

Cr.rm = max{0,2;S|c2p(a-y8|ci|')}, (I5b) 

Ca,n{t) = 2/3^\ci\\c2\, (15c) 
where a, (r,) represents the ith atom (reservoir). 

3. Dynamic analysis 

We apply the asymptotic analysis to atom and reservoir entanglement. Suppose, initially, 
only two atoms are entangled: Co,„2(0) = 2a/3 and C,-i,-2(0) = 0. As ? ^ oo, due to the atom- 
reservoir interaction for each subsystem, atom entanglement inevitably vanishes Caiaiit — ^ 
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Figure 1. Time evolution of concurrence Cf,,,,, (solid curve) and C,.,r, (dashed curve) in non- 
Markovian regime, with the initial state being (fT3l l. for the cases of AjW =0.2 (up plots), (a) 
ff = 1/ V2, (b) a = VTO/5, (c) a = 1/4 and AjW =0.1 (bottom plots), (d) a = 2V2/5, (e) 
a= 1/V5, (f)a= V3/5. 



oo) = and reservoir entanglement necessarily tends to a final value Crir^it ^ oo) = 2a[i. 
This indicates that the information initially stored in the atoms is completely transferred to 
the reservoirs, irrespective of the non-Markovian memory effect. 

In Figure [B we show the time evolution of concurrence between two atoms (solid curve) 
and two reservoirs (dashed curve) as a function of the dimensionless quantity t = At for 
six typical values of the parameter a, namely, a = 1/V2, VTO/5, 1/4 (A/W =0.2) and 
2^/5, 1/V5, V3/5 (^/W^=0.1). 

In Figure \T\ (a), reservoir entanglement performs huge damped oscillations before it 
reaches its final value with no ESD or ESB phenomena and this holds when a > /3. In 
(b), atom entanglement can revive from its last death [fT?l and reservoirs can suddenly be 
entangled in a period of time . In (c), similar to the Markovian case IfTOl . atom entanglement 
vanishes permanently after a finite time, while reservoir entanglement presents sudden birth 
and death even without the revival of atom entanglement, which is totally different from the 
Markovian case where reservoir entanglement just increases monotonically up to a stationary 
value. According to (b) and (c), the revival of atom entanglement does not necessarily indicate 
the sudden death of reservoir entanglement, and vice versa. 

Besides these distinctive non-Markovian behaviors, we point out that entanglement of 
atom pairs can revive before, simultaneously or even after the disentanglement of reservoirs, 
as in Figure [U (d), (e) and (f). According to the expressions of Ca^cii and C^rj, (I15al) and 
i\5b\i respectively, the coincidence of revival and death means 2yS|ci|^(Q' -yS|c2p) > and 
2yS|c2p(Q' - /3\ci\^) < where the two equalities hold at the same time and 2a = /3. We 
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Figure 2. Time evolution of concurrence for different partitions in strong non-Markovian 
limit: Ca^a2 (long-short-dashed curve), C,.j,-, (long-short-short-dashed curve), Ca^r^ (long- 
dashed curve) and Co,,-, (solid curve), with the initial state being (fTsT l and a - Ij VTO, 
AjW - 0.1. The short-dashed curve depicts the population |ci(f)|^. 



can easily prove that revival ahead of (after) death requires 2a > jS {la < /3). These 
counterintuitive phenomena do not rely on the explicit expressions of jcip or [ciP and thus 
are independent of the relative strength of atom-reservoir interaction, that is, the ratio A/W. 

Taking the symmetry of the composite system into account, we analyze the bipartite 
entanglement of the partitions: a^ ® ai, r^ ® r2, ai (8 ri and ai (8 ri in strong non-Markovian 
limit (A/W = OA), as shown in Figure[2l In the initial period of time, due to the atom-reservoir 
couplings, Ca^ai diminishes and Cain arises. According to (I15cl) . Ca,r, has a maximal value fi^ 
with |ci(OI = 1/ V2 and if 2a < a^® aj and ai <S> r2 would already have been disentangled 
at that time. With the successive decay of atoms, the reservoir entanglement emerges and 
evolves to its maximal value la/S with the two-reservoir state being a |0),-j lO)^^ +/5 |l)ri |l)r2- 
In Markovian case, the transmission of information comes to an end. 

However, in non-Markovian case, the information transferred to the reservoirs is fed back 
to the atoms due to the memory effect, which means the excited atomic state will arise again. 
Assuming 5 |ci (01^ /dt = 0, we get sm(dt/2 + 0) = (valley) or tan(<i?/2 + 6) = d/A (peak) 
with 6 = tan~^{d/A). For the peaks, |ci(?„)P = exp(-2nnA/d) (n e N) and = 2nn/d. We 
analyze the behavior of expression 2/3 |cip (a -/3 |c2p) between the peak and valley of |ci(OP 
and find that, the revival of atom entanglement is governed by the inequality 

1 -exp(-2nnA/d) 

a > =. (16) 

^Jl + [1 -exp{-2nnA/d)]^ 

For a given atom-reservoir interaction A/W, the occurrence of revival means 

1 - exp(-27r^/J) 

o: > I (IV) 

^Jl + [1 -exp{-2nA/d)]^ 

If a violates the criterion, Ca,a2 f^il^ to exhibit revival as in Figure [I](c). Furthermore, with a 
given initial condition a , the exact number of revivals is 

na = [{d/2nA)\n(/3/(/3-a))], (18) 
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where y = [x] is the Gaussian function which represents the maximal integer smaller than or 
equal to x. The interval of revivals can be approximately estimated by 

tr ^ tM - ti = In Id, (19) 

which is independent of the initial condition a and largely determined by the transition 
strength W. Therefore, if a satisfies (flTl) . the period of revivals is the same. Here, we note that 
when a > \l V2, the atom entanglement performs damped oscillations and n can be infinite 
as in Figure [H (a). 

In addition, it is not difficult to derive that the sudden death of reservoir entanglement is 
regulated by 

expi-2nnA/d) 

a < (20) 
yl + exp(-4nnA/d) 

Similar to the atom entanglement, the occurrence of BSD means 

exp(-27r^/J) 

a < , (21) 

y 1 + Qx^{-AnAld) 

If a is big enough, there is no BSD between the two reservoirs as in Figure [H (a) and (b). The 
number of BSD is 

nr = [{dl2nA)UI3la)l (22) 

where [■] is the Gaussian function. Based on (fTSl ) and (l22l) . we can see that when a < /3 < 2a, 
revivals of atom entanglement are more than deaths of reservoir entanglement and 2a = 
equal, 2a < /3 less. Though this is visible from the equations, it is not so obvious from our 
physical intuition. 

Back to our analysis of bipartite entanglement, if the memory effect is strong enough, 
that is, the ratio A/W is very small, |ci(?i)| = exp{-nA/d) is bigger than 1/ V2 and C„|,-| will 
rise again to its maximal value. If a satisffes (fTTl) and (|2TI) . the reservoir entanglement will be 
disentangled and the atom reservoir entanglement will revive, as in Figure [2l This means the 
information is transferred from two reservoirs to each subsystem and then to two atoms. This 
depletion-feedback process continues for some time with a damping of amplitudes because 
the memory effect is finite and the atoms will decay inevitably. 

Bntanglement for other partitions is shown in Figure [31 Curve (IV) depicts the 
multipartite entanglement between the four effective qubits, which is defined by multipartite 
concurrence Cn [23]. Other bipartite entanglement is obtained through /-concurrence [2M, 
as curves (I)-(III), (V) and (VI), which coincides with the concurrence in the pure two-qubit 
case. We compile the partitions initially entangled in Figure [3] (a) and disentangled in (b). 
Here, we note that the partition (ai (8) rO (g) (^2 <S) r2) has constant entanglement, specifically 
2ap, despite of the memory effect and information transmission, which serves as a benchmark 
of our entanglement analysis. has the same value at f = and t ^ oo, showing complete 
entanglement transfer from atom pairs to reservoirs, as the asymptotic analysis indicates. 
Although long time evolution of entanglement has the same tendency with the Markovian 
case [IQJ, transient entanglement for different partitions endures huge damped oscillations 
due to the feedback of information. 



Non-Markovian Dynamics of Entanglement for Multipartite Systems 

(a) (b) 





Figure 3. Time evolution of entanglement for different partitions in non-Markovian regime: 
(I) (ai(g)ri)(8)(fl2®'"2), (n)fli®(ri (gia2®'"2), (HI) (ai ® r2)®(fl2®n), (IV) ai®ri (gifl2®'"2, (V) 
(fli ® 02) ® ('"1 ® '"2), (VI) ri (g) (fli (g) 02 ® r2), with the initial state being ( fTsT l and a - Ij VTO, 
AjW = 0.2. 



Our simulation conditions A/W = 0.2 and 0.1 can be well realized within the current 
experimental level f25\ . A more intuitive interpretation of these results will be given in detail 
in the following section. 

4. Explanations via quasimode Hamiltonian 

Memory effect plays an important role in our analysis of non-Markovian dynamics. To 
separate it from the damping effect, we convert the true mode Hamiltonian with Lorentzian 
spectral density into the quasimode form. Applying the method in fil2ll . the quasimode 
Hamiltonian of the single-body system can be given by 



^0 ^memoiy ^damping 5 (23) 



with 



\J — ( 



Ho = 0JQO-+O-- + ojQO^a + Ac"'"(A)c(A)dA, (24a) 



^^memory = W{a+a + a.a^), (24b) 

(a^c(A) + ac+(A))dA, (24c) 

CO 

where c^(A), c(A) are the creation and annihilation operators of the continuum quasimode with 
frequency A and other parameters are the same as before. The particular conversion relations 
are given in the Appendix in [ETl . Analogous to the procedure in section II, the evolution of 
the single-body system is 

\m) = Ca{t) \e) |0)„, 10),, + c,„{t) \g)\l),n 10),. (25) 
+ c,-(0 1^) |0),J1),,, 
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Figure 4. Time evolution of population in strong non-Markovian regime for quasimode 
Hamiltonian: |Cfl(f)p (solid curve), \c,„(t)\^ (short-dashed curve) and |Cr(f)P (long-dashed 
curve), with the initial state being |e) |0),„ |0)r' and A/W = 0.1. 



with Ca(t) = Ci(0, Cm(t) = b{t)exp(-ia)ot) and cXt) = (1 - |Ca(OP - |c,„(OI^)'''^- This means the 
pseudomode in section II is just the discrete quasimode and the quasimodes (both discrete and 
continuum) are our previous reservoir. We plot the time evolution of the population in Figure 
|4]and explain how the memory effect and damping effect induce |ci(OP to perform damped 
oscillations. 

According to (|23]) - (|24cl) . the atom only interacts with the discrete mode and their 
coupling coefficient is just the transition strength W; whereas, the discrete mode interacts 
with a set of continuum modes and their coupling strength contains only the constant width 
of Lorentzian spectral density A. If we let the atom and the discrete mode be a new system 
and the continuum modes be the external environment, the behavior of the new system will 
be exactly Markovian. Therefore, transmission of the exciton from atom to reservoir is a 
two-step process: first, a photon is created in a discrete (cavity) mode via the atom-discrete 
mode interaction; second, this photon is annihilated and a photon is created in a continuum 
(external) mode via the discrete-continuum mode coupling. The memory effect stems from 
the finite life span of the photo in the cavity [r2l[T3]|. Thus, the reabsorbing phenomenon that 
causes the oscillatory entanglement, as shown in Figure [H and [2l only exists in the first step. 
In FigureHl the alternative peaks and valleys of |Ca(OP and |c,„(OP demonstrate that the energy 
of exciton exchanges between the two states \e) \0),„ \0),-, and |^) |1)„, |0)^ periodically. 

Assuming d\Cm(t)\^ /dt = 0, we get s'm(dt/2) = (the valley) or tan(dt/2) = d/A 
(the peak). The valley of |Cm(OP coincides with the peak of Ic^XOP at t = 2nn/d when 
d |Cr(OP /dt = 0, which makes |c,-(Op act like a staircase curve. This means that if the exciton 
only exists in the atom, the damping process does not happen. The damping effect happens 
only when the photon escapes into the continuum modes and never comes back. In our case, 
there is no mechanism, such as the dipole-dipole interaction [211, to protect the exciton from 
escape. So the damping effect is unavoidable and the memory effect is just a finite-time 
phenomenon, which indicates that atom entanglement can never be completely reconstructed 
and the information will be totally transferred to the reservoirs, which fits the results in Figure 



□ and [31 
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5. Conclusions 

We study the dynamical evolution of a couple of two-level atoms interacting with independent 
structured reservoirs. We find that reservoir entanglement exhibits sudden birth, death and 
revival phenomena. The revival of atom entanglement does not necessarily indicate the 
disentanglement of reservoirs, and vice versa. Applying the quantitative analysis, we derive 
the criteria for the revival (death) phenomena and prove that the atom entanglement can 
revive before, simultaneously or after the sudden death of reservoir entanglement, which is 
independent of the relative strength of atom-reservoir interaction. Besides, by studying the 
bipartite and multipartite entanglement for different partitions, we present a comprehensive 
interpretation of the information transmission within this composite system in non-Markovian 
regime. Our results and conclusions are desirable in the implementation of various optical 
schemes [|26l l27l l28l |29ll for the preparation and manipulation of entanglement. 
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